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(57) Abstract: A contact sensitive 
device (10) comprising a member 
(12) capable of supporting bending 
waves, a first sensor (16) mounted on 
the member for measuring bending 
wave vibration in the member, 
the sensor (16) determining a first 
measured bending wave signal 
and a processor which calculates 
information relating to a contact on 
the member (12) from the measured 
bending wave signal, the processor 
applying a correction based on the 
dispersion relation of the material 
of the member (12) supporting the 
bending waves, characterised in that 
the device (10) comprises a second 
sensor (16) to determine a second 
measured bending wave signal which 
is measured simultaneously with the first measured bending wave signal and the processor calculates a dispersion corrected function 
of the two measured bending wave signals which is selected from the group consisting of a dispersion corrected correlation function, 
a dispersion corrected convolution function, a dispersion corrected coherence function and other phase equivalent functions to 
determine information relating to the contact. 
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TITLE: CONTACT SENSITIVE DEVICE 

10 

DESCRIPTION 

15 TECHNICAL FIELD 

The invention relates to contact sensitive devices. 
BACKGROUND ART 

Visual displays often include some form of touch 
sensitive screen. This is becoming more common with the 

20 emergence of the next generation of portable multimedia 
devices such as palm top computers. The most established 
technology using waves to detect contact is Surface Acoustic 
Wave (SAW) , which generates high frequency waves on the 
surface of a glass screen, and their attenuation by the 

25 contact of a finger is used to detect the touch location. This 
technique is n time-of -flight" , where the time for the 
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disturbance to reach one or more sensors is used to detect the 

location. Such an approach is possible when the medium behaves 

in a non-dispersive manner i.e. the velocity of the waves does 

not vary significantly over the frequency range of interest. 

5 In contrast in WO01/48684 to the present applicant, a 

contact sensitive device and method of using the same are 

proposed. The device comprises a member capable of supporting 

. bending wave vibration and a sensor mounted on the member for 

measuring bending wave vibration in the member and for 

10 transmitting a signal to a processor whereby information 
relating to a contact made on a surface on the member is 
calculated from the change in bending wave vibration in the 
member created by the contact. 

By bending wave vibration it is meant an excitation, for 

15 example by the contact, which imparts some out of plane 
displacement to the member. Many materials bend, some with 
pure bending with a perfect square root dispersion relation 
and some with a mixture of pure and shear bending. The 
dispersion relation describes the dependence of the in-plane 

20 velocity of the waves on the frequency of the waves. 

Two types of contact sensitive device are proposed, 
namely a passive sensor in which bending wave vibration in the 
member is only excited by the contact and an active sensor in 
which the contact sensitive device further comprises an 

25 emitting transducer for exciting bending wave vibration in the 
member to probe for information relating to the contact. In 
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the active , sensor, information relating to the contact is 

calculated by comparing the response of waves generated by the 

emitting transducer in the absence of a contact to the 

response caused by the mechanical constraint of the presence 

5 of a contact. 

Bending waves provide advantages, such as increased 

robustness and reduced sensitivity to surface scratches, etc. 

However, bending waves are dispersive i.e. the bending wave 

velocity, and hence the "time of flight", is dependent on 

10 frequency. In general, an impulse contains a broad range of 
component frequencies and thus if the impulse travels a short 
distance, high frequency components will arrive first. This 
effect must be corrected. 

In WO01/48684, a correction to convert the measured 

15 bending wave signal to a propagation signal from a non- 
dispersive wave source may be applied so that techniques used 
in the fields of radar and sonar may be applied to detect the 
location of the contact. The application of the correction is 
illustrated in Figures la to Id. 

20 Figure la shows an impulse in an ideal medium with a 

square root dispersion relation and demonstrates that a 
dispersive medium does not preserve the waveshape of an 
impulse. The outgoing wave (60) is evident at time t=0 and the 
echo signal (62) is spread out over time, which makes a 

25 determination of an exact contact position problematic. 
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In a non-dispersive medium such as air, a periodic 

variation of the frequency response is characteristic of a 

reflection, and is often referred to as comb filtering. 

Physically, the periodic variation in the frequency response 

5 derives from the number of wavelengths that fit between the 

source and the reflector. As the frequency is increased and 

the number of wavelengths fitting in this space increases, the 

interference of the reflected wave with the outgoing wave 

oscillates between constructive and destructive. 

10 Calculating the Fourier transform of the dispersive 

impulse response of Figure la produces the frequency response 
shown in Figure lb. The frequency response is non-periodic 
and the periodic variation with wavelength translates to a 
variation in frequency that gets slower with increasing 

15 frequency. This is a consequence of the square root 
dispersion in which the wavelength is proportional to the 
square root of the inverse of frequency. The effect of the 
panel on the frequency response is therefore to stretch the 
response as a function of frequency according to the panel 

20 dispersion. Consequently, a correction for the panel 
dispersion may be applied by applying the inverse stretch in 
the frequency domain, thus restoring the periodicity present 
in the non-dispersive case. 

By warping the frequency axis with the inverse of the 

25 panel dispersion, Figure lb may be transformed into the 
frequency response for the non- dispersive case (Figure 1c) in 
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which the frequency of excitation is proportional to the 

inverse of the wavelength. This simple relationship translates 

the periodic variation with decreasing wavelength to a 

periodic variation with increasing frequency as shown in 

5 Figure lc. 

Applying the inverse Past Fourier Transform (fft) to the 
trace of Figure lc produces an impulse response shown in 
Figure Id which is corrected for dispersion and where the 
clear reflection is restored. As is shown in Figure Id any 

10 particular waveshape of an impulse is preserved in time since 
the waves travelling in a non- dispersive medium have a 
constant velocity of travel, independent of their frequency. 
Accordingly, the task of echo location is relatively straight 
forward. The outgoing wave (50) is evident at time t=0, 

15 together with a clear reflection (52) at 4ms. The reflection 
(52) has a magnitude which is approximately one-quarter of the 
magnitude of the outgoing wave (50) . 

The procedure described is not applicable if the impulse 
has occurred at an unknown time t Q and the distance x from the 

20 response to an initial impulse may only be calculated if the 
impulse occurs at t Q = 0 . 

It is an object of the present invention to provide an 
alternative contact sensitive device which uses bending wave 
vibration for extracting information relating to the contact. 
25 DISCLOSURE OF INVENTION 

According to one aspect of the invention, there is 
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provided a contact sensitive device comprising a member 
capable of supporting bending waves, a first sensor mounted on 
the member for measuring bending wave vibration in the member, 
the first sensor determining a first measured bending wave 
5 signal and a processor which calculates information relating 
to a contact on the member f rom . the measured bending wave 
signal, the processor applying a correction based on the 
dispersion relation of the material of the member supporting 
the bending waves, characterised in that the device comprises 

10 a second sensor to determine a second measured bending wave 
signal which is measured simultaneously with the first 
measured bending wave signal and the processor calculates a 
dispersion corrected function of the two measured bending wave 
signals which is selected from the group consisting of a 

15 dispersion corrected correlation function, a dispersion 
corrected convolution function, a dispersion corrected 
coherence function and other phase equivalent functions to 
determine information relating to the contact. 

According to a second aspect of the invention, there is 

2 0 provided a method of determining information relating to a 
contact on a contact sensitive device comprising the steps of 
providing a member capable of supporting bending waves and a 
first sensor mounted on the member for measuring bending wave 
vibration in the member, determining, using the sensor, a 

25 ■ first measured bending wave signal characterised by providing 
a second sensor mounted on the member to determine a second 
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measured bending wave signal, measuring the second measured 

bending wave signal simultaneously with the first measured 

bending wave signal, calculating a dispersion corrected 

function of the two measured bending wave signals which is 

5 selected from the group consisting of a dispersion corrected 

correlation function, a dispersion corrected convolution 

function, a dispersion corrected coherence function and other 

phase equivalent functions and processing the measured bending 

wave signals to calculate information relating to the contact 

10 by applying the dispersion corrected function. 

The following features may be applied to both the device 

and the method with the processor being adapted to provide 

many of the calculations or. processing steps of the method. 

The dispersion corrected function may be calculated as 

15 follows: 

calculate W x {coi) and W 2 (ct))* which are the Fourier 
transformation and complex conjugate Fourier transformation of 
the two measured bending wave signals W x (t) and W 2 {t) ; t 
represents time co is 2nf where f is frequency. 
20 calculate a first intermediate function W x (<a)W*(co) ; 

calculate a second intermediate function M(a>) which is a 
function of W x (®)W 2 \o)) ; 

apply a frequency stretching operation f{co), as 
described above in relation to WO01/48684, to M(co) to give the 
25 dispersion corrected correlation function: 
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G{t) = — ^M[f((o)]&vp(iQ)t)dco . 

The intermediate function M(co) may simply be W x {cd)W^{co) 
which gives a standard dispersion corrected correlation 
function. Alternatively, M(co) may be a function which modifies 
5 the amplitude but not the phase of ^(gj)# 2 *(o)) to give a phase 
equivalent function to the standard dispersion corrected 
correlation function. Since the phase equivalent function and 
the standard dispersion corrected correlation function have 
the same phase properties, they have a maximum at the same 
10 position. The phase information in the measured bending wave 
signals may be used to acquire information about the contact 
in particular the location thereof. The location may be 
calculated from the time at which the maximum in the functions 
occurs . 

15 M(co) may be selected from the following functions which 

all yield phase equivalent functions to the standard 
dispersion corrected correlation function: 

a) Ujp^pS^M. 

|i^(<y)r;(c»)| 

Thus M(co) may normalise the amplitudes of W^oa)Wl((o) to 
20 unity to yield a normalised dispersion corrected correlation 
function otherwise known as a dispersion corrected coherence 
function. 



8 



WO 03/005292 PCT/GB02/03073 

Thus M{co) may act on the amplitudes of W^cojW^CD) to 

yield a dispersion corrected correlation function with a 
modified peak shape. 

5 c) Jl^Co^) = ^C^^C^f'I^C^)'^^)!] where <p(x) is a real valued 

function 

Thus M(a>) may apply a general modification to yield a 
phase equivalent function having a different amplitude to the 
standard correlation function. 

10 d) M(o)) = W x (<D)W2{coi)i//((D) where i/r(co) is a real valued 

function 

Thus M(a>) may apply a general frequency- dependent scaling 
to yield a phase, equivalent function having a different 
amplitude to the standard correlation function. Such a 
15 scaling is also known as emphasis. 

Alternatively, M(g>) may be the function D(co) which is 
the Fourier transformation of the correlation function D(t): 

D(o>) is mathematically equivalent to W x (m)W*(a)) and may 
20 be arrived at without calculating W^co) and W 2 (co)*. This is an 
alternative method to calculating the standard dispersion 
corrected correlation function. The steps are calculate D(t); 
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calculate D(co) and apply a frequency stretching operation to 

arrive at the dispersion corrected correlation function: 

One advantage of using the dispersion corrected 
correlation function is that it is applicable in situations 
where the precise time, t 0 , at which a contact occurred is not 
known. This is because an offset t 0 (i.e. t Q *0) in the response 
functions is represented as an additional factor exp(io)t Q ) in 
the Fourier transformations, ^(^)and W 2 ((d) which cancels in 
the intermediate function W^o^W^ico) . 

A transducer may act as both the first and second sensor 
whereby, the dispersion corrected correlation function is an 
autocorrelation function. The autocorrelation function may be 
calculated applying the same steps for the dispersion 
corrected correlation function using Vf\(f)= W 2 (t) . 

The frequency stretching operation may be determined as 
described in relation to Figures la to Id by Q((o) = w\f(<a)\ in 
which f{a>) is chosen so that Q(a>) represents an imaginary 
response of a medium where bending waves of all frequencies 
travel with identical phase velocity. For example, for bending 
mode vibrations that follow a square -root relation of 

wavevector, k, to angular frequency, k = c4co , the frequency 
stretching operation /(#) = vCVtf> gives an imaginary response, 
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Q(t) , of a medium with constant arbitrary phase velocity v . C 

is a constant and is defined by C=(|j/B) A l/4 where ji=mass per 
unit area, and B=bending stiffness. 

The inverse Fourier transform, Q(t) may be examined to 
5 determine the distance to the site of the contact. Cl(t) , is 
centred at a value t x that is proportional to the distance 
between the site of the contact and each sensor. 
The dispersion corrected correlation function may also be used 
to determine the precise location of the contact. For 

10 example, for the function, G(t) , the centre of the correlation 
function is located at t = (x x -x 2 )/v where x l and x 2 are the 
distances from the contact point to the respective vibration 
sensors. The dispersion corrected correlation function may 
therefore be used to determine a difference in path-length 

15 (i.e. x x -x 2 ) between the contact site and the two sensors. 
This quantity, x x -x 2t defines a hyperbola of possible 
locations of the contact on the panel surface. If a third 
sensor is provided, a second dispersion corrected correlation 
function, from a different combination of sensor positions, 

20 provides a second path-length difference, eg. x l -x z . 
Alternatively, a second pair of sensors may be mounted to the 
member to provide a second dispersion corrected correlation 
function. In many cases (for example when the sensors are 
positioned in the four corners of a rectangular panel) , the 

25 two hyperbolic curves defined by the two path-length 
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differences have a unique intersection point which determines 

unambiguously the location of the contact . 

Another advantage of using the dispersion corrected 

correlation function is in the treatment of waves reflected 

5 from boundaries of the member. The reflected waves create the 

effect of virtual sources which are located outside the 

physical boundaries, specifically, at the contact site 

reflected in the axes of the boundaries. The impulse generated 

by a contact may show features corresponding to 'echoes' which 

10 arrive after the direct waves of the impulse. By applying the 
dispersion correction to an impulse at f = 0, discrete 
reflections may be resolved as peaks in the corrected impulse 
response. The location of these peaks may be used as 
additional information for determining the location of the 

15 contact. This may be particularly effective if the 
sensitivity, or acceptance, of the sensors is not independent 
of the direction of the incoming wave since the measured 
signal corresponding to the direct path from the contact may 
be weaker than the signal due to reflected waves . 

2 0 The dispersion- corrected autocorrelation function may be 

applied to locate the path-length difference between the 
sensor and real and virtual sources of bending waves arising 
from reflections at the boundaries. Such information may help 
the determination of the contact location. 

25 Ref lected" waves may be suppressed by placing an absorbing 

material in contact with the edges of the member. 

12 
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Alternatively, the processor may be adapted to remove the 

contribution of reflected waves from the measured bending wave 

signal. This may be achieved on the basis that in a dispersive 

medium, i.e. one with a dispersion relation of the form 

5 k^C-fco , low-frequency components travelling along the direct 
path may arrive after the first high-frequency reflected 
waves. This processor may be adapted to provide a moving 
average low-pass filtering operation on the original measured 
bending wave signal where the width of the averaging window 

10 varies locally with the time coordinate as Atcct 2 . 

The measured bending wave signal may be transformed to a 
response in which the signals due to the direct and reflected 
wave both oscillate with constant but different periods. For 
example, a measured bending wave signal W(t) , from a contact 

15 at time f = 0 may be transformed on the time axis by using 

In the function U(t) the signal from a sharp 
impulse propagating on a dispersive member oscillates with a 
constant period. Furthermore, the oscillations from the direct 
wave have a longer period than those of the reflected waves. 

20 The contribution of reflected waves can therefore be removed 
by a low pass filter operation on U(t) . The response may then 
be transformed back to linear units of time, as W'(t) = U(l/t) . 
The procedure may be generalised to other forms of the 
dispersion relation. 

25 For a perfect rejection of reflected waves, it is 
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necessary to know the dispersion relation in the member, the 

time t 0 at which the contact occurred, and the distance x 
between the contact site and the sensor. However, generally, 
only the first is known. Thus the processor may be adapted to 
5 provide estimates or substitutes for t 0 and x which may be 
substituted into a calculation to remove the reflected waves. 
For example, an estimate for t 0 may be the time at which the 
contact was first detected by any sensor on the member, i.e. 
the time at which the measured signal first passes a 

10 predetermined threshold. The distance x may be set as the 
distance between the relevant sensor and the furthest point on 
the member or the maximum dimension (e.g. diameter) of the 
member. Use of these substitutes should leave the direct-path 
signal intact; Estimates of t 0 and x obtained by other means 

15 may also be used. 

There may be multiple (i.e. n) sensors on the member and 
thus the number of distinct correlation functions is n(/z-l)/2 . 
The processor may be adapted to create a mapping function 
which maps the surface of the member for each correlation 

20 function whereby the dispersion corrected correlation 
function, G(t), is swept over the surface such that all co- 
ordinates of a given path-length difference, Ax, take the 
value G(Ax/v) . The product of the entire set of mapped 
correlation functions may then be calculated and the location 

25 of the contact may be determined from the maximum co-ordinate. 
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This method has the desired property of increasingly rejecting 

spurious contributions from reflections as the number of 

sensors is increased. 

Each sensor and its associated circuitry may have 

5 identical or different phase properties. Differences in the 

phase properties may be caused by mechanical inconsistencies 

in the manufacture and mounting of the sensors, or by finite 

tolerances of component values in the associated amplification 

and filtering circuits. If the phase properties of each 

10 sensor differ, the dispersion corrected correlation function 
may be convolved with some asymmetric function. This 
asymmetric function may lead to incorrect measurements of 
path-difference difference since the position of peaks in the 
dispersion corrected correlation function is no longer 

15 independent of the phase properties of the sensor. 

The processor may thus be adapted to perform the 
following steps: 

a) estimate a convolution correction coefficient & 12 (o)) from 



j*; 2 (*>)= 



2 0 where {^(0)} and ^T*2,y(o)} are the Fourier transformation and 
complex conjugate Fourier transformation of two measured 
bending wave signals §V hj {i)\ and \jV 2J (t)) and {Axj\ is the path- 
length difference; 

b) calculate the dispersion corrected correlation function 
25 with phase correction from: 

15 
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The phase correction has the effect of deconvolving the 
dispersion corrected correlation function such that the 
position of the maximum is consistent with the position of the 
contact . 

5 The processor may further be adapted to include in the 

determination procedure any available information about where 
the contact can be expected. This may be particularly useful 
in situations where the position of a contact may not be 
unambiguously determined by the dispersion corrected 

10 correlation functions, e.g. when reflected waves interfere 
with the direct -wave signature or less than three sensors are 
used. For example, if the member is an input device for a 
graphical user interface where the user is presented with a 
choice of x buttons' to press, it may be useful to assume that 

15 any contact on the member occurs within the discrete areas 
corresponding to the buttons. 

Alternatively, a map of the probability at which a 
contact is likely to occur and which is based on the expected 
behaviour of the user may be used. The device may comprise a 

20 software application with a graphical user interface (GUI) 
which interacts with the operating system by means of an 
application program interface (API) in which the API is 
adapted to generate the probability map. The probability map 
may be based on the location, size, and frequency of use of 
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objects presented by the graphical user interface. The 

probability map may also be based on information about the 

relative likelihood of the various GUI elements being 

activated. 

5 The information in the mapped correlation functions may 

be combined with the probability map to give higher 
reliability. The expected probability may also be a further 
input to a neural -net which accepts information derived from 
the sensors. 

10 The member may comprise a raised pattern on its surface 

whereby a contact drawn across the surface provides a variable 
force to the member to generate bending waves in the member. 
The pattern may be periodic, or quasi -periodic with a 
statistically well-defined spatial distribution of 

15 undulations. The processor may be adapted to determine the 
contact position by using knowledge of the periodicity of the 
pattern and the fact that the interval between impulses 
represents the time in which a contact, which may be provided 
by a stylus, has travelled to an adjacent feature of the 

20 pattern. The pattern may be random whereby a contact 
travelling over the surface of the member generates a random 
bending wave signal. 

The use of random surface pattern may be used 
independently of the dispersion correction correlation 

25 function. Thus according to another aspect of the invention, 
there is provided a contact sensitive device comprising a 
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member capable of supporting bending waves, a sensor mounted 

on the member for measuring bending wave vibration in the 

member to determine a first measured bending wave signal and a 

processor which calculates information relating to the contact 

from the measured bending wave signal from the sensor, 

characterised in that a surface of the members comprises a 

raised pattern whereby a contact drawn across the surface 

provides a variable force to the member to generate bending 

waves in the member. 

The device may comprise sensing means to determine a 
second measured bending wave signal which is measured 
simultaneously with the first measured bending wave signal and 
the processor may calculate information relating to the 
contact from a dispersion corrected correlation function of 
the two measured bending wave signals. The dispersion 
corrected correlation function is described above and thus 
features of the first and second embodiments may be applied to 
this embodiment. The dispersion corrected correlation 

function is particularly useful if the measured bending wave 
signals have a large bandwidth centred on a frequency which 
gives a phase -velocity of bending waves in the member which is 
much greater than the maximum lateral velocity of the contact. 

According to another aspect of the invention, there is 
provided a contact sensitive device comprising a member 
capable of supporting bending waves, a sensor mounted on the 
member for measuring bending wave vibration in the member to 
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determine a first measured bending wave signal and a processor 

which calculates information relating to a contact from the 

measured bending wave signal from the sensor, characterised in 

that the device comprises at least a second sensor to 

5 determine a second measured bending wave signal which is 

measured simultaneously with the first measured bending wave 

signal and the processor optimises a product of a set of 

corrected impulse response measurements from each sensor to 

determine information related to the contact. 

10 The information calculated may be the time t 0 at which 

contact occurs and thus the form of the dispersion relation, 

k(oo) , and the distance, x, separating each sensor from the 

contact site must be known. The corrected impulse response 

measurement may be calculated by using the following steps: 

15 1) Calculate the Fourier transform W(co) of a measured bending 
wave signal W(t) ; 

2) Calculate an equivalent response, W'(t) , from a notional 
sensor positioned at the contact site by 

WXco) = W(a>)exp[ik(a))x\ . 

3) Calculate the inverse Fourier transform of W\co) to provide 

20 function W'(t) . 

The product is thus YLffi^) in which the function W'if) 
shows an initial impulse from the contact which is without 
dispersion and which is centred at t 0 , the time of the 
impulse. Features in W(t) due to edge reflections will appear 
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later in W'(t) but will not be corrected for dispersion in the 

same way as for the direct impulse which is thus more easy to 
identify. 

The product shows a strong peak due to the direct 
impulse, and a weak background noise. By taking the product we 
reinforce the signal coincident at t Q in all W}(t) , whereas the 
information due to reflections is uncorrelated and suppressed. 
This procedure is therefore a reliable way of determining f 0 . 

Conversely, the information relating to the contact may 
be the location of the contact, i.e. the co-ordinate r and 
thus the sensor-contact distances, Xj and the time t Q for 

which the maximum value of YLffi^o) is obtained must be known. 

This optimisation process may involve iterative 
refinement of estimates for r and t 0 . The initial estimate may 
be derived from impulse response functions whose high temporal 
frequency components have been suppressed, with the result 
that the estimate is easy to obtain (fewer local maxima) , but 
of low spatial precision. Further iterations may introduce 
progressively higher frequency components as the estimate is 
refined. 

The invention therefore provides two complementary 
methods of determining the contact position: the method of 
dispersion-corrected correlation functions, and the method of 
maximising J^J RP} f (f 0 ) . 

The following characteristics may apply to all 

20 
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embodiments of the invention. The device may comprise means 

for recording measured bending wave signals from the or each 

sensor over time as the contact moves across the member. The 

measured bending wave signals may be recorded as time-series 

5 data, i.e. a sequence of values measured at different times. 

The time -series data may be analysed as a sequence of 

short sections or 'frames' of data, which may have widths or 

durations of 10ms. The processor may apply a window function 

to the frames of data. Window functions are well known in the 

10 art - see for example Press W.H. et al . , Numerical Recipes in 
C, 2 Ed., Cambridge University Press 1992, Ch. 13.4. The 
processor may be adapted to extract information on the contact 
which has been averaged over the duration of the frame, e.g. 
the mean positional co-ordinate of the moving contact. The 

15 processor may apply the dispersion corrected correlation 
technique to calculate the mean path length difference for 
each frame of data to give a sequence of path- length 
differences over time. 

Alternatively, the time-series data may be analysed using 

2 0 adaptive filters such as those described in Grant PM et al 
"Analogue and Digital Signal Processing and Coding" , Ch 10 
(1989) . The adaptive filter may correct for a convolution 
function which relates the time-series data from the first and 
second measured bending wave signals. The convolution function 

25 changes slowly over time as the contact moves and is dependent 
on the position of the contact. 
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The adaptive filter may calculate a convolution 

correction which when applied to the measured bending wave 

signals produces a signal as close as possible to the 

measurement input. The first measured bending wave signals 

5 form a first time -series data which may be sent to the 

reference signal of the adaptive filter and the second 

measured bending wave signals form a second time- series data 

which may be delayed, preferably by the maximum expected width 

of the correlation function, before being sent to the primary 

10 input of the adaptive filter whereby the convolution 
correction may be calculated. The processor may apply the 
dispersion correction to the convolution function to give a 
function whose maximum is determined by the difference between 
the path-length from contact to the sensor and the path-length 

15 from the contact to the sensing means. 

The information calculated may be the location of the 
contact or may be other information, e.g. pressure or size of 
the contact. The information relating to the contact may be 
calculated in a central processor. The sensors may be mounted 

2 0 at or spaced from an edge of the member. The sensors may be 
in the form of sensing transducers which may convert bending 
wave vibration into an analogue input signal. 

The member may be in the form of a plate or panel . The 
member may be transparent or alternatively non- transparent, 

25 for example having a printed pattern. The member may have 
uniform thickness. Alternatively, the member may have a more 
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complex shape, for example a curved surface and/or variable 

thickness. The processor may be adapted for complex shaped 

members by providing an adaptive algorithm such as a neural 

net to decipher the contact location from the bending wave 

5 signal received by the sensor. 

The device may be a purely passive sensor with the 

bending wave vibration and hence the measured bending wave 

signals being generated by an initial impact or by frictional 

movement of the contact. Alternatively, the device may be an 

10 active sensor and thus the device may comprise an emitting 
transducer. The transducer may have dual functionality, namely 
acting as an emitting transducer and a sensor. The member may 
also be an acoustic radiator and bending wave vibration in the 
member may be used to generate an acoustic output. 

15 Measurements of the bending wave signal due to a contact 

may be contaminated by the bending waves due to the audio 
. signal, particularly when the audio signal is similar to the 
bending wave signal generated by the contact. The effect may 
be minimised by ensuring the frequency band of the audio 

20 signal differs from and does not overlap the frequency band of 
the measurements from the sensors and sensing means. The 
audio and measured signals may thus be filtered, for example, 
the audio band may be limited to frequencies below 20kHz, and 
the vibration measurements may be limited to frequencies above 

25 20kHz. 

The device may be a dual active and passive sensor and 
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may be adapted to switch between active and passive sensing 
modes depending on whether contact is applied to the device. 
The device may cycle between resting in passive sensing mode 
when no contact is detected, switching to active mode sensing 
5 when a contact is applied and returning to passive sensing 
mode once the contact is removed to wait for further contacts. 
This may be advantageous to avoid the power drain when the 
device is in active mode. 

The contact may be in the form of a touch from a stylus 

10 which may be in the form of a hand-held pen. The movement of a 
stylus on the member may generate a continuous signal which is 
affected by the location, pressure and speed of the stylus on 
the member. The stylus may have a flexible tip, e.g. of 
rubber, which generates bending waves in the member by 

15 applying a variable force thereto. The variable force may be 
provided by tip which alternatively adheres to or slips across 
a surface of the member. As the tip moves across of the 
member a tensile force may be created which at a certain 
threshold, causes any adhesion between the tip and the member 

20 to break, thus allowing the tip to slip across the surface. 

Alternatively, the contact may be in the form of a touch 
from a finger which may generate bending waves in the member 
which may be detected by passive and/or active sensing. The 
bending waves may have frequency components in the ultrasonic 

25 region (>20 kHz) . Passive sensing is therefore sensitive to 
contacts with both fingers and styli. 
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When the device is acting as an active sensor, i.e. with 
an emitting transducer generating an excitation signal, the 
contact may exert a non- linear force on the member so as to 
generate harmonics of the excitation signal. The processor may 
5 comprise signal processing devices to isolate the excitation 
signal from the harmonics so that the harmonics may used to 
determine the contact position in a similar maimer to passive 
sensing. The harmonics effectively constitute a source of 
bending waves from the contact site. 

10 The or each emitting transducer or sensor may be a bender 

transducer which is bonded directly to the member, for example 
a piezoelectric transducer. Alternatively, the or each 
emitting transducer or sensor may be an inertial transducer 
which is coupled to the member at a single point. The inertial 

15 transducer may be either electrodynamic or piezoelectric. It 
may be possible to use audio transducers which are already in 
place as sensing and/or emitting transducers. 

A contact sensitive device according to the invention may 
be included in a mobile phone, a laptop or a personal data 

20 assistant. For example, the keypad conventionally fitted to a 
mobile phone may be replaced by a continuous moulding which is 
touch sensitive according to the present invention. In a 
laptop, the touchpad which functions as a mouse controller may 
be replaced by a continuous moulding which is a contact 

25 sensitive device according to the invention. The moulding may 
be implemented as a mouse controller or other alternatives, 
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e.g. a keyboard. Alternatively, the contact sensitive device 
may be a display screen, e.g. a liquid crystal display screen 
comprising liquid crystals which may be used to excite or 
sense bending waves. The display screen may present 
5 information relating to the contact. 
BRIEF DESCRIPTION OF DRAWINGS 

The invention is diagrammatically illustrated, by way of 
example, in the accompanying drawings, in which: - 

Figures la to Id are a graphic illustration of a method 
10 of dispersion correction according to the prior art, in which 
Figure la is a graph of a dispersive impulse response showing 
response in arbitrary units against time. Figure lb is a graph 
of a dispersive frequency response showing response in 
arbitrary units against frequency. Figure lc is a graph of a 
15 non- dispersive frequency response showing response in 
arbitrary units against frequency. Figure Id is a graph of a 
non-dispersive impulse response showing response in arbitrary 
units against time; 

Figure 2 is a plan view of a contact sensitive device 
20 according to the present invention; 

Figure 3 is a perspective view of a first device 
incorporating passive touch sensing; 

Figure 4 is a perspective view of a second device 
incorporating passive touch sensing; 
25 Figure 5 is a block diagram of a processing algorithm for 

the passive sensing of Figures 3 and 4; 
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Figure 6 is a perspective view of first device 

incorporating active touch sensing; 

Figure 7 is a perspective view of a second device 

incorporating active touch sensing; 

Figure 8 is a flow chart showing a method of calculating 

the path- length difference between the contact and two 

measurement points using the dispersion corrected correlation 

function; 

Figure 8a is a schematic plan view of a device to which 
the method of Figure 8 is applied; 

Figure 8b is a graph of dispersion corrected correlation 
function against time; 

Figure 9 is a flow chart showing a first method of 
removing reflections from a measured bending wave signal, 

Figure 10 is a flow chart showing a second method of 
removing reflections from a measured bending wave signal. 

Figure 11a and lib are schematic perspective and plan 
views of a touch sensitive device according to another aspect 
of the invention; 

Figure 12 is a schematic block diagram of a processing 
algorithm which may be used for the device of Figure 11. 

Figure 12a is a graph showing the combined transfer 
function H(f) against frequency (f) for the filter and 
amplifier of Figure 12; 

Figure 12b is a graph of a measured bending wave signal 
against time; 
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Figure 13 is a flow chart of the steps for obtaining an 
empirical phase correction; 

Figure 13a is a plan view of a grid defined on a member 
for use in the method of Figure 13; 
5 Figure 14 is a touch sensitive device according to 

another aspect of the invention; 

Figure 15 is a circuit diagram of an adaptive noise 
canceller which may be used in the various devices; 

Figure 16a is a schematic block diagram of a contact 
10 sensitive device which also operates as a loudspeaker; 

Figure 16b is a method of separating audio signal and 
measured bending wave signal in the device of Figure 16a; 

Figure 17 is a flow chart showing a method of calculating 
the contact location using the dispersion corrected auto- 
15 correlation function; 

Figure 17a is a schematic plan view of a device to which 
the method of Figure 17 is applied; 

Figure 17b is a graph of dispersion corrected auto- 
correlation function against time, and 
20 Figure 18 is a block diagram showing how an adaptive 

filter may be used to calculate information relating to the 
contact . 

Figure 2 shows a contact sensitive device (10) comprising 
a transparent touch sensitive plate (12) mounted in front of a 
25 display device (14) . The display device (14) may be in the 
form of a television, a computer screen or other visual 
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display device. A stylus (18) in the form of a pen is used 
for writing text (20) or other matter on the touch sensitive 
plate (12) , 

The transparent touch sensitive plate (12) is a member, 
e.g. an acoustic device, capable of supporting bending wave 
vibration. Three transducers (16) are mounted on the plate 
(12) . At least two of the transducers (16) act as sensors or 
sensing means and are thus sensitive to and monitor bending 
wave vibration in the plate. The third transducer (16) may 
also be a sensing transducer so that the system corresponds to 
the passive contact sensitive device of Figure 3 or Figure 4 . 

Alternatively, the third transducer may be an emitting 
transducer for exciting bending wave vibration in the plate so 
that the system corresponds to the active sensor of Figure 5. 
In the Figure 6 or Figure 7 embodiment, the active sensor may 
act as a combined loudspeaker and contact sensitive device. 

Figures 3 and 4 are more detailed illustration of two 
contact sensitive devices (32,33). The contact sensitive 
devices (32,33) comprises a member in the form of a panel (24) 
capable of supporting bending wave vibration and three sensors 
in the form of sensing transducers (26) for sensing bending 
wave vibration at their respective mounting points. The 
vibration pattern (28) is created when pressure is applied at 
a contact point (30) . The devices may be considered to be 
passive contact sensitive devices since the devices do not 
comprise an emitting transducer. Thus the bending wave panel 
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vibration in the panel is generated solely by the contact. 

In a passive sensor an impulse in the body of the panel 
(24) starts a bending wave travelling towards the edge of the 
panel (24) . The bending wave is detected by the three sensing 
5 transducers (2 6) mounted eguidistantly around the edges as in 
Figure 3 or by the three sensing transducer mounted on a 
surface of the panel (24) but spaced from the edges of the 
panel (24) as in Figure 4. The measured bending wave signals 
are processed to determine the spatial origin and force 
10 profile of the applied impulse. 

Figure 5 shows a possible implementation for the 
processing of the bending wave information sensed at each 
sensing transducer (26) of Figure 3 or Figure 4. .In Figure 5, 
the bending waves in the panel are sensed by three sensing 
15 transducers (26) . The sensing transducers (26) measure 
analogue bending wave signals Wi(t), W 2 (t) and W 3 (t) which are 
transmitted to a multiplexing analogue to digital converter 
(ADC) (54) . The resultant digital input signal is transmitted 
to the central processor (34) from which information (58) 
2 0 relating to the location and profile of the contact impulse is 
determined. 

Figures 6 and 7 are more detailed illustrations of 
alternative combined touch sensitive and audio devices 
(35,37) . The devices each comprise a panel (24) capable of 
25 supporting bending wave vibration and an emitting transducer 
(31) for exciting bending wave vibration in the panel (24) . 
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The device (35) in Figure 6 further comprises two sensing 
transducers (26) for sensing bending wave vibration at their 
respective mounting points whereas the device (37) in Figure 7 
comprises only one sensing transducer (26) . The vibration 
5 pattern (28) is interrupted when pressure is applied at a 
contact point (30) . The devices may be considered to be active 
contact sensitive devices since the devices comprise an 
emitting transducer (31) . 

In Figure 6, the sensing and emitting transducers (26,31) 

10 are spaced equidistantly around the edges of the panel (24) 
whereas in Figure 7, the sensing and emitting transducers 
(26,31) are distanced from the edges of the panel (24) and are 
mounted to a surface thereof. 

Figure 8a shows an embodiment having two sensors (102) 

15 mounted on a member (100) to which a contact is applied at a 
contact location (104) . Figure 8 shows a method of calculating 
the dispersion corrected correlation function to reveal the 
difference in path length between the contact location (104) 
and the sensors (102) . The method comprises the following 

20 steps: 

(a) Measure two bending wave signals W x (t) and W 2 (t) ; 

(b) Remove reflections from the measured signals and calculate 
W({t) and Wi(t),e.g. by using the method set out in Figure 9; 

(c) Calculate the Fourier transform of W x '(t) and W{{t) to 
25 arrive at W x (a>) and W 2 (co) and hence the intermediate function 
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W x {a)) Wl((D) ; where W^(co) is the complex conjugate Fourier 
transform. 

(d) and (e) at the same time as performing steps (a) to (c) , 
the frequency ' stretching operation f{(o)-vC4o is calculated 
5 using the predetermined panel dispersion relation k = c4& . 

(f) W x {cS) and W 2 {oj) and f{co)-vC4c^ are combined to arrive at 
the dispersion corrected correlation function: 

G{f)~ fVil/W^I/^Jexp^O^ ; and 
27V J - fl0 

(g) the dispersion corrected correlation function is plotted 
10 against time with a peak occurring at time t 12 as shown in 

Figure 8b; 

(h) AX12 is calculated from t i2 ; Ax i2 is the path-length 
difference between the path lengths x x and x 2 from the first 
and second sensors to the contact . 

15 (i) Axi 2 is used to calculate the location of the contact. 

Alternatively at step (e) , the dispersion corrected 
correlation function with phase correction <j>i 2 set out below 
may be used. The calculation of <j> 12 is explained in Figure 13. 

20 Figure 9 shows a method of removing reflections from an 

impulse measurement where the impulse occurs at t=0. The 
method comprises the following steps: 
i) Measure a bending wave signal W x (i) ; 
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ii) Transform the signal on the time axis by using C/(r)=^(l/r) ; 

iii) A low pass filter is applied, for example, U' (t) as shown 
to remove all reflected signals. The -constant C from the 
predetermined panel dispersion relation k-C-[a> is used to 

5 define the width of the convolution function At; 

iv) The response may then be transformed back to linear units 
of time, as W'(t) = U(l/t) . 

Figure 10 shows an alternative method for removing 
reflections from an impulse measurement where the impulse 
10 occurs at t=0. The method comprises the following steps: 

i) Measure a bending wave signal W x (t) ; 

ii) Estimate the distance x between each sensor and the 
location of the contact; 

iii) Use the estimate x and a predetermined dispersion 

15 relation k = Cjco to define an averaging window N(t, t 1 ) 

iv) Apply the averaging window N(t, t') to the bending wave 
signal W x (f) to remove the effect of reflections. 
N(t,t') is an example of an averaging window with a 

Gaussian shape. A rectangular window can be implemented to 
20 give computational efficiency. Multiple applications of a 

rectangular window may be desirable; a large number of 

repeated applications of a rectangular window will produce a 

similar result to N(t,t'). 

Figures 11a and lib show a contact sensitive device (80) 
25 comprising a rectangular member (82) capable of supporting 
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bending waves and four sensors (84) for measuring bending wave 
vibration in the member. The sensors (84) are in the form of 
piezoelectric vibration sensors and are mounted on the 
underside of the member (82), one at each corner. A foam 
5 mounting (86) is attached to the underside of the member and 
extends substantially around the periphery of the member. The 
foam mounting (86) has adhesive surfaces whereby the member 
may be securely attached to any surface. The foam mounting may 
reduce the reflections from the edge of the member. 

10 Two sets of path length differences, Axi 2 =Xi-x 2 and 

Ax 3 4=x 3 -X4 are calculated as described in Figure 8; xi is the 
distance from each sensor to the contact. As shown in Figure 
lib the hyperbolae (85) defined by the two path- length 
differences are plotted and the location of the contact is the 

15 intersection (87) of two hyperbolae. 

Figure 12 is a schematic diagram illustrating the 
implementation of a processing algorithm in the device of 
Figure 11. The sensors (84) measure analogue bending wave 
signals W x (t) # W 2 (t) and W 3 (t) which are passed through an 

20 amplifier and anti-aliasing (low-pass) filter (88) . The 
amplitude of the combined transfer function H(t) of the 
amplifier and anti-aliasing filter is shown in Figure 12a. The 
filtered signals are converted into digital signals by a 
digitiser (90) and stored in a first-in-first-out buffer 

25 having finite length. The buffer comprises two stores, a pre- 
trigger and a post-trigger store (92,93) for signals measured 
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before and after the detection process is triggered 
respectively. 

The central processor (94) determines information relating 
to the location and profile of a contact of the member by the 
following steps: 

a) The central processor performs a threshold test which is 
illustrated in Figure 12b. The measured- bending wave 
signal (96) is compared to a predetermined threshold value 
(98) . When the measured signal passes the threshold 
value, the detection process is triggered. 

b) An array of time-series digital input signals is 
transferred from the buffer to the processor. The signals 
include measurements taken before and after the detection 
process is triggered so that a digital measure of the 
entire waveform of the impulse from the contact is 
reconstructed. 

c) The processor shifts the waveform according to the 
estimate of t 0 so that to is set to zero. 

d) The processor removes the effect of reflections from the 
digitised signal as described above with the estimate of 
t 0 as zero and the estimate of x taken as the diagonal 
length of the member. 

e) The processor applies further processing, in particular 
calculating the dispersion corrected correlation function 
for each diagonally opposed pair of sensors and 
calculating information relating to the contact. 
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The further processing applied by the processor may also 
comprise applying a pre -determined phase correction to the 
dispersion corrected correlation function. This may be 
calculated as shown in Figure 13 . 

a) Define a grid - for example, for the embodiment shown in 
Figures 11a and lib which has a member (82) which, four 
sensors (84) , this may be done by defining points (71) by 
the vector {rj}. In this example the grid has 8 rows and 7 
columns so the grid is defined by a set of 56 vectors [ri, 

b) A user taps on the first point (71) in the grid defined by 
ri and the first and second bending wave signals {W i/ j(t)} 
and {W 2 ,j (t) } are measured by the first and second sensor of 
each pair of sensors; 

c) Step (b) is repeated until the user has tapped on each 
point in the grid; 

d) Calculate the Fourier transforms of the bending wave 
signals; 

e) Calculate the empirical phase correction: 



where Axj is the difference between the path lengths, x if j & x 2f j 
from the first and second sensors to the contact. The path 
lengths are known from the grid coordinates. 

Figure 14 shows a contact sensitive device (70) 
comprising a member (72) capable of supporting bending waves 
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and three sensors (64) mounted on the member for measuring 

bending wave vibration in the member. A surface of the member 

(72) comprises a raised pattern (66) which is a periodic 

pattern of raised crossed lines. A stylus (78) is drawn across 

the surface along a path (74) and as it crosses a line of the 

pattern it generates bending waves (76) in the member. 

Figure 15 shows an adaptive noise canceller for example 
as described in "Widrow et al Adaptive Noise Cancelling: 
Principles and Applications, Proceedings of the IEEE, Vol 63 
No 12 pp 1692 (1975)". The adaptive noise canceller comprises 
an adaptive filter (40) which takes the initial audio signal 
as the reference input. The canceller may be used to remove 
the contribution of the audio signal from the output of a 
vibration sensor before any further processing occurs. The 
adaptive filter shown is one example of an adaptive filter 
that can be applied to this task. 

Figure 16a shows a contact sensitive device which also 
operates as a loudspeaker. Figure 16b shows a method for 
partitioning the audio signal and measured signal into two 
distinct frequency bands so that the contribution of the audio 
signal to the processed measured signal is suppressed. The 
device comprises a member (106) in which bending waves are 
generated by an emitting transducer or actuator (108) and the 
contact . The emitting transducer applies an audio signal to 
the member (106) to generate an acoustic output. Before being 
applied to the member, the audio signal is filtered by a low 
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pass filter (112) which, as shown in Figure 16b, removes the 
audio signal above a threshold frequency f 0 . 

As shown in Figure 16b, the contact generates a signal 
which has a power output which is substantially constant over 
5 a large frequency band. The signal from the contact and the 
audio signal sum to give a combined signal which is passed 
through a high pass filter (114) to remove the signal above 
the threshold frequency f 0 . The filtered signal is then 
passed to a digitiser (116) and onto a processor (118) . 

10 Figure 17a shows an embodiment having a single sensor 

(120) mounted on a member (100) to which a contact is applied 
at a contact location (104) . Bending waves are reflected from 
the edge of the member and create an image of a virtual source 
which is at location (122) . Figure 17 shows a method of 

15 calculating the dispersion corrected auto -correlation function 
to reveal the contact location (104) . The method comprises 
the following steps: 

(a) Measure one bending wave signal W x (t) ; 

(b) Calculate the Fourier transform of W((t) to arrive at 



(c) at the same time as performing steps (a) and (b) , the 



20 W x {a>)i 



frequency stretching operation /(o)) = vCVfi) 



is calculated using 



the predetermined panel dispersion relation k = c4co . 



(d) W x {a>) and f{co)^vC4w 



are combined to arrive at the 



25 dispersion corrected auto -correlation function: 
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(g) the dispersion corrected correlation function is plotted 
against time with peaks occurring at time t n and -tn as shown 
in Figure 17b; 

5 (h) Ax xl is calculated from tn; Ax i2 is the path-length 
difference between the path lengths x x and xi« from the first 
and second sensors to the contact. 

(i) AX12 is used to calculate the location of the contact. 

Figure 18 shows the adaptive filter of Figure 15 may also 
10 be used to calculate the location of the contact from a device 
comprising two sensors. In general, adaptive filters contain a 
finite-impulse-response (FIR) filter. A FIR filter is 
equivalent to a convolution operation, with some convolution 
function 0 (t) . 

15 The signals W x (t) measured by the first sensor are sent 

to the adaptive filter (40) and the signals W 2 (t) measured by 
the second sensor are sent to a delay unit (41) . The delay 
unit delays the signals from the second sensor, preferably by 
the maximum expected width of the convolution function. The 

20 delayed signals are then sent to the primary input of the 
adaptive filter. The adaptive filter continually updates the 
convolution function so that an estimate, W 2 (t) , of the primary 
signal input, W 2 (t) , may be obtained from the reference input, 
Wi(t) . The convolution operation is defined as follows: 
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The location of the contact is calculated in the 
processor by the following steps: 

a) Extract the FIR convolution function from the internal 
5 memory of the adaptive filter. 

b) Calculate the Fourier transform of the FIR convolution 
function. 

c) Apply the frequency stretching operation /(a>) = vCV# ; 

d) Calculate the inverse Fourier transform to arrive at F(t) . 
10 F(t) is a phase equivalent of the dispersion correction 

correlation function G(t) and thus the Fourier transforms of 
G(t) and F(t) have equal phase but not necessarily the same 
amplitude. Accordingly, the location of any peaks in the time 
domain for F(t) and G(t) are the same and thus the location of 
15 the contact may be calculated from the peak of F(t) as 
described above for G(t) . 
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CLAIMS 

1. A contact sensitive device comprising a member capable of 
supporting bending waves, a first sensor mounted on the member 
for measuring bending wave vibration in the member, the sensor 
5 determining a first measured bending wave signal and a 
processor which calculates information relating to a contact 
on the member from the measured bending wave signal, the 
processor applying a correction based on the dispersion 
relation of the material of the member supporting the bending 

10 waves, characterised in that the device comprises a second 
sensor to determine a second measured bending wave signal 
which is measured simultaneously with the first measured 
bending wave signal and the processor calculates a dispersion 
corrected function of the two measured bending wave signals. 

15 which is selected from the group consisting of a dispersion 
corrected correlation function, a dispersion corrected 
convolution function, a dispersion corrected coherence 
function and other phase equivalent functions to determine 
information relating to the contact. 

20 2. A contact sensitive device according to claim 1, wherein 
the first sensor acts as the second sensor whereby the 
calculated dispersion corrected function is an autocorrelation 
function. 

3. A contact sensitive device according to claim 1 or claim 
25 2, comprising a second pair of sensors to determine two 
additional measured bending wave signals from which a second 



41 



WO 03/005292 PCT/GB02/03073 
dispersion corrected function is calculated. 

4. A contact sensitive device according to claim 3, wherein 
the processor determines from each dispersion corrected 
function a first difference in path- length between the contact 

5 site and each of the first and second sensors and a second 
difference in path-length between the contact site and each of 
the second pair of sensors and determines the location of the 
contact from the first and second differences in path-length. 

5 . A contact sensitive device according to any preceding 
10 claim, comprising absorbers mounted around at least part of 

the periphery of the member to absorb reflected waves. 

6. A contact sensitive device according to any preceding 
claim, wherein the processor is configured to remove the 
contribution of reflected waves from the measured bending wave 

15 signal. 

7. A contact sensitive device according to claim 6, wherein 
the processor comprises a low-pass filtering operator which 
operates on the measured bending wave signal and which 
comprises an averaging window having a width which varies 

20 locally with time. 

8 . A contact sensitive device according to claim 6 or claim 
7, wherein the processor is configured to provide an estimate 
for the distance between the contact site and each sensor, the 
estimate being substituted into a calculation to remove the 

25 reflected waves. 

9. A contact sensitive device according to any preceding 
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claim, comprising multiple sensors on the member whereby- 
multiple dispersion corrected functions are determined. 

10. A contact sensitive device according to claim 9, wherein 
the processor is configured to create a mapping function which 
maps the surface of the member for each dispersion corrected 
function. 

11. A contact sensitive device according to any one of the 
preceding claims, wherein the member is an acoustic radiator 
and an emitting transducer is mounted to the member to excite 
bending wave vibration in the member to generate an acoustic 
output . 

12. A contact sensitive device according to claim 11, 
comprising means ensuring that the acoustic output and 
measured bending wave signals are in discrete frequency bands. 

13. A contact sensitive device according to claim 12, 
comprising an adaptive noise canceller for removing the 
contribution of the acoustic output from the measured bending 
wave signal . 

14. A contact sensitive device according to any one of the 
preceding claims, wherein the member is transparent. 

15. A contact sensitive device according to any one of the 
preceding claims, wherein the processor estimates a 
convolution, correction coefficient which is applied to the 
dispersion corrected function thereby compensating for phase 
differences between the sensors. 

16. A contact sensitive device according to any one of the 
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preceding claims, wherein the member comprises a raised 

pattern on its surface whereby a contact drawn across the 

surface provides a variable force to the member to generate 

bending waves in the member. 

5 17. A contact sensitive device according to claim 16, wherein 

the pattern is random whereby random bending wave vibration is 

generated in the member. 

18. A contact sensitive device according to any one of the 
preceding claims, comprising means for recording sets of first 

10 and second measured bending wave signals from the or each 
sensor over time as the contact moves across the member. 

19. A contact sensitive device according to claim 18, wherein 
the processor is adapted to analyse the measured bending wave 
signals as a sequence of frames of data. 

15 20. A contact sensitive device according to claim 18, wherein 
the processor is adapted to extract information on the contact 
which has been averaged over the duration of the frame. 

21. A contact sensitive device according to any one of the 
preceding claims, wherein the processor comprises an adaptive 

2 0 filter which calculates a convolution function between the set 
of first measured bending wave signals and the set of second 
measured bending wave signals. 

22. A contact sensitive device according to claim 21, wherein 
the processor is adapted to use a dispersion corrected 

25 convolution function to calculate information about the 
contact . 
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23. A method of determining information relating to a contact 
on a contact sensitive device comprising the steps of 

providing a member capable of supporting bending waves 

and a first sensor mounted on the member for measuring bending 
5 wave vibration in the member, 

determining, using the sensor, a first measured bending 

wave signal characterised by- 
providing a second sensor mounted on the member to 

determine a second measured bending wave signal, 
10 measuring the second measured bending wave signal 

simultaneously with the first measured bending wave signal, 

calculating a dispersion corrected function of the two 

measured bending wave signals which is selected from the group 

consisting of a dispersion corrected correlation function, a 
15 dispersion corrected convolution function, a dispersion 

corrected coherence function and other phase equivalent 

functions and 

processing the measured bending wave signals to calculate 
information relating to the contact by applying the dispersion 
20 corrected function. 

24. A method according to claim 23, wherein the dispersion 
corrected correlation function is calculated by the following 
steps: 

calculate the Fourier transformation and complex 
25 conjugate Fourier transformation of the two measured bending 
wave signals; 
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calculate an intermediate function by multiplication of 

the Fourier transformation and complex conjugate Fourier 

transformation of the two bending wave signals; 

combine a frequency stretching operation and subsequent 

inverse Fourier transformation with the intermediate function 

to give the dispersion corrected correlation function. 

25. A method according to claim 23, wherein a phase 
equivalent function is calculated by the following steps: 

calculate the Fourier transformation and complex 
conjugate Fourier transformation of the two measured bending 
wave signals; 

calculate a first intermediate function by multiplication 
of the Fourier transformation and complex conjugate Fourier 
transformation of the two bending wave signals; 

calculate a second intermediate function which is a 
function of the intermediate function; 

combine a frequency stretching operation and subsequent 
inverse Fourier transform with the second intermediate 
function to give the phase equivalent function. 

26. A method according to claim 25, wherein the second 
intermediate function is calculated by normalising the first 
intermediate function to give the dispersion corrected 
coherence function. 

27. A method according to any one of claims 24 to 26, 
comprising removing the effect of reflections before 
calculating the or each intermediate function. 
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28, A method according to any one of claims 23 to 27, 

comprising using the dispersion corrected function to 
determine a difference in path-length between the contact site 
and each sensor and hence to determine the location of the 
5 contact . 

29. A method according to any one of claims 23 to 28, 
comprising providing a transducer which acts as both the first 
and second sensor and calculating a dispersion corrected 
autocorrelation function. 

10 30. A method according to any one of claims 23 to 29, 
comprising providing a second pair of sensors to determine two 
additional measured bending wave signals from which a second 
dispersion corrected function is calculated. 

31. A method according to claim 30, comprising determining 
15 from each dispersion corrected function a first difference in 

path-length between the contact site and each of the first and 
second sensors and a second difference in path-length between 
the contact site and each of the second pair of sensors and 
determining the location of the contact from the first and 
20 second differences in path-length. 

32. A method according to any one of claims 23 to 31, 
comprising removing the contribution of reflected waves from 
each measured bending wave signal by operating a low-pass 
filtering operation on the measured bending wave signals. 

25 33. A method according to any one of claims 23 to 31, 
comprising providing multiple sensors on the member whereby 
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multiple dispersion corrected functions are determined. 

34. A method according to claim 33, comprising creating a 
mapping function which maps the surface of the member for each 
dispersion corrected function. 

35. A method according to any one of claims 23 to 34, 
comprising calculating a convolution correction coefficient 
and applying the convolution correction coefficient when 
calculating the dispersion corrected function to compensate 
for any phase differences between the sensors. 

36. A method according to any one of claims 23 to 35, 
comprising recording a set of first measured bending wave 
signals and a set of second measured bending wave signals from 
the or each sensor over time as the contact moves across the 
member . 

37. A method according to claim 36, comprising analysing the 
sets as a sequence of frames of data. 

38. A method according to claim 37, comprising extracting 
information relating to the contact from each frame which has 
been averaged over the duration of the frame. 

39. A method according to claim 37 or claim 38, comprising 
correcting for a convolution function between the sets, 

40. A contact sensitive device comprising a member capable of 
supporting bending waves, a sensor mounted on the member for 
measuring bending wave vibration in the member to determine a 
first measured bending wave signal and a processor which 
calculates information relating to a contact from the measured 
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bending wave signal from the sensor, characterised in that the 
device comprises at least a second sensor to determine a 
second measured bending wave signal which is measured 
simultaneously with the first measured bending wave signal and 
5 the processor optimises a product of a set of corrected 
impulse response measurements from each sensor to determine 
information related to the contact. 

41. A contact sensitive device according to claim 40, wherein 
the corrected impulse response measurement is calculated by 

10 using the following steps: 

Calculate the Fourier transform of the measured bending 
wave signal; 

Calculate an equivalent response from a notional sensor 
positioned at the contact site 
15 Calculate the inverse Fourier transform of the equivalent 

response to provide a function to be optimised. 

42. A contact sensitive device according to claim 41, wherein 
the optimisation process involve iterative refinement of 
estimates of the location of the contact and the time for 

20 which the maximum value of the product is obtained. 

43. A contact sensitive device according to claim 42, wherein 
an initial estimate is derived from impulse response functions 
whose high frequency components have been suppressed. 

44 . A contact sensitive device comprising a member capable of 
25 supporting bending waves, a sensor mounted on the member for 

measuring bending wave vibration in the member to determine a 
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first measured bending wave signal and a processor which 
calculates information relating to the contact from the 
measured bending wave signal from the sensor, characterised in 
that a surface of the members comprises a raised pattern 
5 whereby a contact drawn across the surface provides a variable 
force to the member to generate bending waves in the member. 

45. A contact sensitive device according to claim 44, wherein 
the pattern is periodic, 

46. A contact sensitive device according to claim 45, wherein 
10 the pattern is quasi -periodic with a statistically well- 
defined spatial distribution of undulations. 

47. A contact sensitive device according to any one of claims 
40 to 46, wherein the processor is configured to determine the 
contact position by using knowledge of the periodicity of the 

15 pattern and the fact that the interval between impulses 
represents the time in which a contact has travelled to an 
adjacent feature of the pattern. 

48. A contact sensitive device according to any one of claims 
1 to 22 and 40 to 48, wherein the device is a purely passive 

2 0 sensor with the bending wave vibration and hence the measured 
bending wave signals are generated by an initial impact or by 
frictional movement of the contact. 

49. A contact sensitive device according to any one of claims 
1 to 22 or 40 to 48, wherein the device is an active sensor 

25 comprising an emitting transducer. 

50. A contact sensitive device according to claim 48 or claim 
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49 , wherein the device is a dual active and passive sensor and 

is configured to switch between active and passive sensing 

modes depending on whether contact is applied to the device. 

51. A contact sensitive device according to claim 50, wherein 

the device cycles between resting in passive sensing mode when 

no contact is detected, switching to active mode sensing when 

a contact is applied and returning to passive sensing mode 

once the contact is removed to wait for further contacts . 
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